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[Text]  Abstract 

This  paper  describes  the  design  method  for  an  extended  nozzle  with  gas- 
deployed  skirt  and  double-nested  exit  cone,  and  estimates  the  performance 
gain  of  this  nozzle.  The  results  show  that  the  nozzle  has  significant 
practical  value. 

I.  Introduction 

The  extended  nozzle  is  a  new  nozzle  design  which  has  been  under  development 
over  the  past  20  years  for  use  on  intercontinental  [ballistic]  missiles  and 
the  final  stages  of  space  planes;  it  is  generally  acknowledged  that  this  new 
nozzle  can  significantly  increase  the  range  or  the  payload  of  the  vehicle. 

The  extended  nozzles  in  existence  today  can  be  divided  into  two  types.  One 
is  a  gas-deployed  skirt  design  which  resembles  a  folded  umbrella;  it  is 
actually  a  metallic  skirt  folded  along  the  radial  direction  toward  the  inner 
chamber  of  the  nozzle.  Its  advantages  are;  light  weight,  simple  structure, 
ease  of  deployment,  and  pre-seal  design;  however,  its  application  is  rather 
limited  (the  area  ratio  must  be  greater  than  25) .  The  other  type  is  a 
nested-tube  exit  cone  which  can  be  extended  during  operation;  its  design 
consists  of  a  fixed  exit  cone  (also  called  the  base  nozzle)  and  one  or  two 
extension  nozzles  which  can  be  extended  to  the  designated  operating  position 
by  an  executive  mechanism  (which  may  be  a  rigid  or  a  flexible  actuator  tube) . 
The  nested  type  extended  nozzle  is  structurally  simple,  operationally 
reliable,  and  has  adjustable  expansion  ratio  to  provide  higher  performance 
gain.  But  the  number  of  exit  cones  is  limited  to  a  maximum  of  three. 

Table  1  shows  a  comparison  of  the  performance  of  exit  cone  between  the 
nested  tube  design  and  the  gas-deployed  skirt  design.^  It  is  not  difficult 
to  see  that  one  can  design  a  superior  extended  nozzle  by  incorporating  the 
desirable  features  and  eliminating  the  shortcomings  from  each  of  the  above 
designs. 
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Table  1.  Comparison  of  the  Performance  of  Exit  Cone  Between  the  Nested  Tube  Design  and  the  Gas-Deployed 
Skirt  Design 
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Over  the  past  10  years,  significant  progress  has  been  made  in  this  country  in 
the  development  of  extended  nozzles. 2  Both  the  gas-deployed  skirt  type 
nozzles  and  the  nested-tube-type  nozzles  have  been  tested,  and  a  technology 
base  has  been  established  for  designing  a  composite  type  extended  nozzle. 

II.  Design  of  the  Extended  Nozzle 

1.  Design  Criteria  and  Requirements 

The  design  criteria  and  requirements  for  the  extended  nozzle  include: 
operating  conditions,  engine  performance,  and  weight  and  spatial  constraints. 
They  are  discussed  below: 

1)  Operating  conditions,  i.e.,  the  operating  altitude  and  environmental 
requirements  (such  as  overload,  shock,  vibration,  storate  and  operating 
temperature,  nuclear  radiation  environment); 

2)  Separation  and  operating  mode:  whether  ignition  takes  place  before  or 
after  deployment; 

3)  Weight  and  spatial  constraints;  the  allowable  nozzle  weight  and  the  space 
between  stages  (which  determine  the  retraction  length  and  the  radial  dimen¬ 
sions  of  the  extended  nozzle) ; 

4)  Engine  dimensions,  thrust  level,  combustion-chamber  temperature,  operating 
pressure  and  duration; 

5)  Throat  diameter,  submerge  depth,  exit  half  angle  and  expansion  ratio  of  the 
nozzle; 

6)  Deployment  length,  deployment  time  and  speed; 

7)  Nozzle  swivel  angle,  velocity  and  acceleration  as  well  as  the  swivel 
moment  and  load  position; 

8)  Materials  used  and  their  properties; 

9)  Stability  and  reliability  requirements; 

10)  Economy,  etc. 

2.  Practical  Design  Considerations 

An  extended  nozzle  with  a  gas-deployed  skirt  and  a  double-nested-tube  design 
has  three  segments:  a)  the  base  nozzle  and  the  double-nested  extendable  exit 
cone  (the  latter  has  a  metallic  skirt  which  can  be  folded  in  the  radial 
direction);  b)  the  deplo3nnent  system;  and  c)  the  power-supply  system.  The 
basic  principle  of  this  design  is  to  fully  take  advantage  of  the  results  and 
the  mature  technologies  that  have  been  developed  for  the  gas-deployed  skirt 
and  single/double  rigid-action  swiveling  extended  nozzles.  Figure  1  shows  a 
schematic  diagram  of  the  extended  nozzle  design. 
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Figure  1.  Schematic  Diagram  of  a  Composite  Extended  Nozzle 


Some  of  the  practical  design  issues  will  be  discussed  in  each  of  the  following 
areas; 

(1)  Design  of  the  Base  Nozzle  and  the  Extended  Exit  Cone 


1)  Determination  of  the  number  of  exit  cones.  The  number  of  extended  exit 
cones  is  determined  on  the  basis  of  spatial  constraints,  performance  require¬ 
ments,  and  the  effect  on  thrust  vector  control.  A  double  nested  exit  cone  is 
used  in  this  design. 

2)  Determination  of  the  length  of  the  exit  cone  and  the  exit  diameter.  A 
very  long  base  exit  cone  implies  a  short  extended  exit  cone  whose  effective¬ 
ness  would  be  limited;  on  the  other  hand,  a  very  short  base  exit  cone  will 
result  in  an  excessively  long  and  heavy  extended  exit  cone  and  significantly 
increase  the  weight  of  the  deployment  system.  Therefore,  the  length  and 
diameter  of  the  exit  cone  (including  the  gas-deployed  skirt)  must  be  care¬ 
fully  chosen  to  keep  the  expansion  ratio  relatively  small  and  to  keep  the 
exit  cone  diameter  in  its  retracted  position  to  be  less  than  the  engine 
diameter. 


3)  Design  of  the  gas-deployed  skirt.  The  metallic  skirt  is  either 
mechanically  attached  or  bonded  to  the  second  section  of  the  extended  exit 
cone,  and  folded  along  the  radial  direction.  Clearly,  the  size  of  the  skirt 
is  constrained  by  the  diameter  of  the  exit  cone;  however,  with  careful 
design,  it  is  possible  to  have  the  deployed  exit  diameter  greater  than  the 
engine  diameter,  and  to  keep  the  gas-deployed  skirt  from  extending  until  the 
movable  cone  is  deployed  so  that  the  exhaust  gas  would  not  enter.  Such  a 
design  is  particularly  desirable  on  strategic  missiles. 

k 


4)  Design  of  the  inner  surface.  In  order  to  avoid  particle  collisions,  this 
design  uses  a  special  nozzle  with  an  idealized  surface  for  the  throat  region 
and  for  the  entrance  section  of  the  exit  cone.  To  satisfy  the  spatial  con¬ 
straints  of  the  exposed  part  of  the  nozzle  and  to  meet  the  performance  gain 
requirements,  the  length  of  the  nozzle  is  chosen  to  be  as  short  as  possible 
and  its  expansion  ratio  is  chosen  to  be  as  large  as  possible;  also  the 
shape  of  the  nozzle  surface  is  chosen  to  be  a  third-order  polynomial,  as 
shown  in  Figure  2.  In  its  fully  extended  and  deployed  position,  the  nozzle 
length  is  1010  mm,  and  the  inner  diameter  of  the  nozzle  exit  is  825.1  mm, 
the  expansion  ratio  can  be  as  high  as  45.  For  the  purpose  of  comparison,  the 
figure  also  shows  the  inner  surface  of  a  fixed  exit  cone  with  the  same  exposed 
length. 


X 

Figure  2.  Inner  Surface  of  a  Composite  Extended  Nozzle  in  Its 
Fully  Extended  and  Deployed  Position 

5)  Materials  for  the  outer  shell  and  the  skirt.  The  outer  shell  of  the 
nozzle  is  made  of  aluminum  alloy;  the  gas-deployed  skirt  is  made  of  niobium 
alloy;  in  other  countries,  the  skirt  is  generally  made  of  metallic  actinium 
and  stainless  steel. 3 

6)  Ablation  material  for  the  exit  cone.  Initially,  silicon  oxide  fiberglass 
wrapped  with  cloth  tape  is  used  as  ablation  material;  the  improved  design 
will  use  carbon-carbon  composite  material  to  reduce  the  mass  of  the  exit  cone. 

7)  Sealing  material.  Silicon  rubber  0-rings  are  used  not  only  as  sealing 
material  but  also  as  a  cushion  to  absorb  the  impact  load. 

8)  Locking  mechanism.  A  mechanical  locking  mechanism  is  used,  as  in  the  case 
of  the  single-section  swiveling  extended  nozzle. 

9)  Insulation  material.  In  order  to  prevent  damages  caused  by  gas  leaks,  it 
is  necessary  to  coat  the  actuator  tube,  the  exit  cone,  the  gas-deployed  skirt 
and  the  support  structure  with  a  layer  of  protective  metallic  oxide  (e.g., 
aluminum  oxide  or  silicon  oxide) . 

10)  Buffer  material.  Conventional  sponge-rubber  plates  (the  improved  design 
will  use  soft  graphite)  are  used  as  the  buffer  material  to  reduce  vibration 
and  to  absorb  the  impact  load  during  deployment. 
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11)  Continuity  of  the  inner  surface.  Steps  or  discontinuities  at  the  joints 
of  the  exit  cone  will  cause  the  heat  conduction  coefficient  along  the  wall 
downstream  of  the  steps  to  increase  and  therefore  result  in  severe  ablation 
of  the  nozzle.  In  order  to  obtain  a  better  surface,  these  steps  must  be 
weakened;  however,  based  on  fabrication  and  assembly  considerations,  the  steps 
must  retain  sufficient  thickness  to  avoid  defects  or  damages. 

(2)  Design  of  the  Deployment  System 

A  complete  deployment  system  should  contain  the  following  segments:  the 
positioning  and  directional  guidance  segment,  the  locking  unit,  the  sealing 
unit  and  the  buffer  unit.  The  primary  requirements  for  a  deployment  system 
are:  good  deployment  performance,  accurate  synchronization,  axial  S37mmetry 
and  controllability.  In  designing  the  deployment  system,  one  must  carry  out 
load  calculations  which  take  into  account  the  inertial  load  and  the  pre¬ 
tension  of  seal  rings,  as  well  as  the  thrust  and  aerodynamic  load  after  the 
nozzle  is  deployed.  Aerodynamic  loads  must  also  be  calculated  in  designing 
the  actuator  tube.  Clearly,  the  above  loads  will  not  occur  simultaneously. 

The  Initial  stage  of  deplo3nnent  is  dominated  by  inertial  load  and  aero¬ 
dynamic  load;  the  final  stage  of  deplo5nnent  is  dominated  by  the  pre-tension 
and  the  Increased  thrust  force.  For  each  stage  one  can  determine  the  maximum 
load  for  each  of  the  actuator  tubes. 

This  design  uses  four  nested  actuator  tubes  (a  detailed  description  is  pre¬ 
sented  by  the  designer  in  a  separate  paper) .  Gas  is  passed  through  four 
equal-length  manifolds  into  the  actuator  tubes;  when  the  pressure  Inside  the 
manifold  reaches  a  certain  value,  the  actuator  tube  becomes  pressurized,  and 
begins  to  extend  and  rotate.  In  the  deployed  configuration,  the  travel  of 
the  first  section  of  the  exit  cone  is  100  mm,  and  that  of  the  second  section 
is  200  mm.  Each  actuator  tube  is  equipped  with  a  pulley  system  to  ensure 
synchronized  travel.  At  the  end  of  the  deplo3rment  sequence,  the  mouth  of  the 
exit  cone  is  pressed  against  the  cushion  and  the  seal  ring,  and  it  is  locked 
in  the  deployed  position  by  a  mechanical  locking  device.  After  the  nozzle  is 
deployed  and  before  engine  cut-off,  the  pressure  and  rigidity  inside  the 
actuator  tube  are  maintained  by  the  pressure  in  the  manifold. 

The  metal  skirt  which  is  folded  along  the  radial  direction  is  deployed 
directly  by  the  low-pressure  gas  from  engine  ignition,  and  is  stabilized  in 
the  open  position  by  the  internal  aerodynamic  pressure  of  the  jet  stream  from 
the  exhaust  nozzle. 

(3)  Selection  of  Power  Supply 

During  cold  tests  and  ground  tests,  the  rigid  actuator  tube  is  operated  by 
compressed  air;  however,  when  conditions  permit,  it  can  also  be  powered  by 
the  gas  from  a  gas  generator.  To  develop  this  new  power  supply,  effort 
should  be  initiated  to  study  the  design,  fabrication  and  testing  of  gas 
generators . 
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III.  Estimation  of  Performance  Gain 


The  performance  gain  that  is  achievable  by  the  extended  nozzle  can  be  esti¬ 
mated  as  follows.  First,  the  performance  is  calculated  based  on  the 
expansion  ratio;  this  performance  estimate  is  modified  by  the  additional 
propellant  due  to  reduction  in  the  submerge  depth  of  the  nozzle  and  the 
performance  degradation  caused  by  the  Increase  in  structural  weight.  If  the 
shockwave  loss  in  the  supersonic  flow  caused  by  discontinuities  of  the  inner 
surface  is  neglected,  then  the  net  performance  gain  of  the  extended  nozzle 
over  a  fixed  nozzle  with  the  same  exposed  length  can  be  determined. 

The  mass  of  the  exit  cone  of  the  nozzle  can  be  estimated  by  the  following 
formula;^  a  more  accurate  estimate  can  be  obtained  by  the  method  of 
integration. 


where  p  is  the  density  of  the  material; 

"S’  is  the  average  wall  thickness; 
a  is  the  average  exit  half  angle; 
dj[  is  the  throat  diameter; 

Sexit  is  the  exit  area  ratio  of  the  extended  exit  cone; 

Sentry  i®  entry  area  ratio  of  the  extended  exit  cone; 

AM  is  the  mass  of  the  deplojnnent  system; 

M  is  the  mass  of  the  exit  cone. 

Table  2  shows  the  estimated  performance  gain  of  the  extended  nozzle. 
Specifically,  the  results  for  three  different  designs  of  extended  nozzles 
and  for  the  fixed  nozzle  (with  the  same  exposed  length  and  diffusion  loss) 
are  compared.  It  can  be  seen  from  the  table  that  design  No.  3,  which  uses 
aluminum  alloy  and  carbon-carbon  composite  material,  can  exceed  10s,  which 
is  4s  higher  than  the  performance  gain  of  the  double-nested-type  extended 
nozzle  (without  the  gas  skirt) .  Such  a  design  can  significantly  enhance 
the  range  or  the  payload  of  the  rocket. 

IV.  Conclusion 

The  development  of  a  double-nested-type  extended  nozzle  with  a  gas-deployed 
skirt  is  a  very  complex  task;  there  are  many  issues  that  require  further 
study  in  order  to  arrive  at  an  improved  and  more  reliable  design.  Clearly, 
by  using  new  light-weight,  high-strength  metal  alloys  and  carbon-carbon 
composite  materials,  the  structural  weight  of  this  composite  extended  nozzle 
can  be  significantly  reduced  and  its  specific  impulse  can  be  further 
increased.  Such  a  design  will  undoubtedly  have  wide  range  of  applications 
on  space  planes  and  rockets . 
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Table  2.  Estimate  of  the  Performance  Gain  of  Composite  Type  Extended  Nozzles 


\  Item 

\  name 

Nozzle  section 

NDesign 

\  \ 
Nozzle 

type  \ 

Nozzle  performance 

_ 1 

Exit  cone 

Deployed  skirt 

Total 

mass 

(kg) 

\ 

Throat 

diam. 

(mm) 

Expan¬ 

sion 

ratio 

Theoreti¬ 

cal 

specific 

impulse 

Material 

Material 

Extended 

nozzle 

1 

123 

45 

292.1 

Aluminum, 

silicon 

oxide 

65.95 

Niobium 

alloy 

8.88 

74.83 

123 

45 

292.1 

Aluminum, 

silicon 

oxide 

65.95 

Niobium 

alloy 

8.88 

74.83 

3 

123 

45 

292.1 

Aluminum, 

carbon- 

carbon 

50.87 

Niobium 

alloy 

8.88 

59.75 

Fixed 

nozzle 

1 

123 

13.94 

274.7 

Aluminum, 

silicon 

oxide 

45 

2 

123 

13.94 

274.7 

Aluminum, 

carbon- 

carbon 

35 

8 


Actuator 

tube 

Gas  generator 

Total 

system 

mass 

(kg) 

Specific  impulse  gain  (s) 

Material 

Mass 

(kg) 

Material 

n 

Increase 

in 

theore¬ 

tical 

specific 

impulse 

Increase 
in  pro¬ 
pellant 

Reduction 
in  specific 
impulse  due 
to  increase 
in  struc¬ 
tural  mass 

Net 

speci¬ 

fic 

impulse 

gain 

Steel 

14.97 

Aluminum, 

packed 

powder 

2.14 

91.94 

17.4 

0.534 

9.95 

7.98 

Steel, 

aluminum 

alloy 

7.92 

Aluminum, 

packed 

powder 

2.14 

84.89 

17.4 

0.534 

8.46 

9.47 

Steel, 

aluminum 

alloy 

7.92 

Aluminum, 

packed 

powder 

2.14 

69.81 

17.4 

0.534 

7.38 

10.55 
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[Text]  Abstract 

This  paper  describes  the  physical  process  of  electrical  discharge  inside  an 
orificed  hollow  cathode.  Based  on  the  design  parameters  of  an  8-cm  mercury- 
ion  thruster  and  the  hollow  cathode,  the  pressure,  particle  density,  and  the 
energy-exchange  mean  free  path  of  the  initial  electrons  inside  the  cathode 
are  calculated.  From  these  results,  the  dimensions  of  the  cathode  emitter 
and  the  top  orifice  are  determined,  and  the  cathode  emission  current  is 
estimated.  The  thermal  equilibrium  parameters  and  heater  parameters  of  the 
cathode  are  calculated  and  compared  with  experimental  results.  The  struc¬ 
tural  design  of  the  cathode  is  discussed  according  to  its  operating  mode  in 
the  thruster;  the  discussion  includes  selection  of  materials  and  plasma 
coating  techniques.  In  addition,  the  results  of  long-term  thermal-simulation 
experiments  are  also  presented. 

I.  Introduction 

In  an  8-cm  electron-bombardment-type  mercury-ion  thruster,  two  self- 
sustaining  arc-discharge  thermal  cathodes  are  used.  One  of  the  cathodes 
whose  top  side  is  located  at  the  center  of  the  base  plane  of  the  ion  source, 
is  called  the  discharge-chamber  cathode.  Its  main  function  is  to  initiate 
ion  discharge  and  to  maintain  constant  discharge  from  the  chamber  during 
normal  thruster  operation;  the  required  discharge  current  is  0.5A.  The  other 
cathode,  located  outside  the  ion-source  housing,  is  called  the  neutralizer 
cathode.  Its  function  is  to  inject  an  equivalent  electron  current  into  the 
plasma  beam  in  order  to  keep  the  propulsion  system  electrically  neutral;  the 
required  discharge  current  is  0.55A. 

In  the  current  designs  of  electron-bombardment-type  mercury-ion  thrusters, 
both  cathodes  are  orificed  hollow  cathodes  with  thermal  electron  emitters, 
as  shown  in  Figure  1;  the  structures  of  the  two  cathodes  are  basically 
identical.  However,  since  the  mercury  discharge  through  the  neutralizer 
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cathode  is  one-tenth  the  amount  of  the  discharge-chamber  cathode,  the  dimen¬ 
sions  of  the  top  orifice  and  the  operating  parameters  of  the  two  cathodes  are 
not  the  same. 


Figure  1.  Structure  of  an  Orificed  Hollow  Cathode 

1.  cathode  top  with  center  hole;  2.  internal  emitter;  3.  contact 
electrode;  4.  thermal  screen;  5.  heater;  6.  AI2O3  coating; 

7.  tungsten  coating;  8.  tantalum  sleeve;  9.  ceramic  tube; 

10 .  tantalum  tube 

During  the  life  of  a  thruster,  the  cathode  is  required  to  fire  several 
thousand  times;  during  each  firing,  it  is  subject  to  severe  thermal  shock 
and  the  bombardment  of  ion  discharge.  Therefore,  the  cathode  is  one  of  the 
components  of  the  thruster  that  is  the  most  fragile  and  whose  performance 
deteriorates  the  fastest.  For  this  reason,  special  care  must  be  taken  in 
selecting  the  cathode  parameters  and  in  its  structural  design.  Articles  on 
the  performance,  manufacturing  techniques  and  thermal-simulation  results  of 
the  discharge-chamber  cathode  and  neutralizer  cathode  used  in  an  8-cm, 

5  mN  [milliNewton]  mercury-ion  thruster  can  be  found  in  the  literature;  this 
paper  will  focus  its  discussion  on  the  cathode  design. 

Because  the  orificed  hollow  cathode  with  Internal  emitter  has  many  advantages 
such  as  low  power  consumption,  durability  and  resistance  against  ion  bombard¬ 
ment,  it  has  been  used  successfully  on  the  ground  as  a  replacement  for  a 
metallic  cathode  in  ion-coating  processes  and  in  neutral-beam  injectors;  the 
new  cathode  has  a  much  longer  life  span  and  leads  to  improved  performance  of 
the  equipment.  It  can  also  be  used  as  an  effective  gas-discharge  electron 
source  in  inert-gas  thrusters,  argon-ion  lasers  and  high-intensity  electron 
beams . 

II.  Design  of  the  Inner  Diameter  of  the  Cathode  Tube  and  the  Diameter  of  the 
Top  Orifice 

In  order  to  achieve  effective  ionization  of  the  mercury  atoms  by  the  initial 
electrons  emitted  inside  the  cathode,  and  to  maintain  a  stable  discharge,  it 
is  necessary  to  design  the  cathode  diameter  to  match  the  energy-exchange  mean 
free  path  of  the  initial  electrons,  Xp^*  Thus,  under  a  given  set  of  condi¬ 
tions,  the  value  of  Xpj-  should  be  first  determined,  and  then  the  diameter  of 
the  cathode  tube  can  be  determined  accordingly. 
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1.  Pressure  Inside  the  Cathode  Tube 

In  general,  different  parts  of  the  cathode  tube  have  unequal  pressures 
because  of  the  motion  of  the  mercury  gas  and  the  presence  of  the  plasma. 
However,  since  the  gas  flow  velocity  near  the  small  top  orifice  is  very  low, 
the  pressure  difference  across  the  cathode  primarily  occurs  in  the  immediate 
vicinity  of  the  top  orifice  and  inside  the  orifice  channel.  Therefore,  one 
can  regard  the  pressure  to  be  approximately  constant  inside  the  cathode;  it 
is  referred  to  as  the  static  pressure.  This  pressure  is  a  function  of  the 
mercury  discharge,  the  magnitude  of  the  discharge  current,  and  the  diameter 
of  the  top  orifice;  the  functional  relationship  can  be  represented  by  the 
following  empirical  formula:^ 

Po=J|^(13.7  +  7.82//.)x  133.3x  10'®  (1) 

d  0 

where  Pq  is  the  static  pressure  inside  the  cathode; 

m  is  the  mercury  discharge  (in  equivalent  Amperes) ; 
do  is  the  diameter  of  the  top  orifice; 

Id  is  the  discharge  current  of  the  cathode. 

2.  Density  of  Neutral  Particles  Inside  the  Cathode 

The  total  pressure  inside  the  cathode  is  equal  to  the  sum  of  partial 
pressures  of  the  individual  particles 

P  =  n^kT^  +  nekT^+nikTt  (2) 

where  n  and  T  are  the  density  and  temperature  of  the  particles,  and  the  sub¬ 
scripts  0,  e,  i  refer  to  the  neutral  particles,  electrons  and  ions 
respectively;  k  is  Boltzman's  constant. 

As  pointed  out  earlier,  one  can  make  the  approximations  P  =  Po  and  To  =  T^  = 
Tg,  where  Tg  is  the  cathode  wall  temperature.  For  a  quasi-neutral  plasma, 
i^i  “  ’^e*  equation  (2)  can  then  be  rewritten  in  the  following  form; 


_  Po — Tlekijr e -\-T t) 

kT. 


(3) 


Test  results  show  that  the  value  of  Tg  remains  relatively  constant  over  a 
wide  range  of  values  of  Ip,  A  and  do,  and  is  equal  to  0.71  ±  0.1  eV  (8236K),1 
and  the  value  of  ng  is  only  a  fraction  of  1  percent  of  no.  Thus,  even  if 
Tg  =  8236K,  Tg  =  1273K,  the  second  term  of  equation  (3)  is  only  5  percent 
of  no*  As  a  first-order  approximation,  equation  (3)  can  be  written  as; 


(4) 
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3.  Energy-Exchange  Mean  Free  Path  of  the  Initial  Electrons 

Inside  the  cathode,  there  is  a  high  frequency  of  occurrence  of  non-elastic 
collisions  between  the  initial  electrons  and  mercury  atoms  and  elastic 
collisions  between  the  initial  electrons  and  Maxwell  electrons.  The  energy 
emitted  into  the  plasma  by  the  Initial  electrons  is  consumed  primarily  to 
excite  and  to  ionize  the  mercury  atoms;  the  energy  exchange  process  involves 
converting  kinetic  energy  into  internal  energy.  The  elastic  collisions 
between  the  initial  electrons  and  Maxwell  electrons  Involve  only  exchange 
of  kinetic  energies,  and  have  an  effect  on  the  energy-exchange  mean  free 
path  of  the  initial  electrons.  Elastic  collisions  also  take  place  between 
the  initial  electrons  and  mercury  atoms  and  ions,  but  the  frequency  of 
occurrence  is  several  orders  of  magnitude  lower  than  that  of  the  above- 
mentioned  collisions;  therefore,  their  effect  can  be  ignored.  The  effective 
energy-exchange  mean  free  path  of  the  initial  electrons  can  be  expressed  as: 


where  is  the  mean  free  path  of  non-elastic  collisions  between  the  initial 
electrons  and  mercury  atoms; 

Xe  is  the  mean  free  path  of  elastic  collisions  between  the  initial 
electrons  and  Maxwell  electrons. 

The  initial  electrons  emitted  by  the  internal  emitter  are  accelerated  by  the 
plasma  sheath  potential  along  a  radial  trajectory.  Inside  the  cathode  tube, 
they  oscillate  in  the  radial  direction  (because  they  cannot  overcome  the 
repulsive  field  to  return  to  the  emitter)  until  they  collide  with  other 
particles;  at  the  same  time,  they  also  migrate  downstream  under  the  action  of 
the  weak  axial  electric  field,  and  exit  the  tube  through  the  top  orifice.  As 
a  result,  most  of  the  energy-exchange  activities  of  the  particles  occur  in 
the  vicinity  of  the  orifice;  this  is  the  region  where  plasma  density  is  the 
highest,  and  the  bombardment  of  ions  on  the  emitter  is  the  most  intense, 
hence  the  temperature  is  also  the  highest.  Test  results  show^"^  that  the 
distribution  of  plasma  density  decreases  exponentially  along  the  distance 
upstream  from  the  orifice;  the  temperature  of  the  emitter  reaches  a  peak 
value  within  a  narrow  section  at  the  downstream  end  of  the  tube,  and 
decreases  linearly  with  distance  upstream.  When  the  pressure  inside  the 
cathode  exceeds  400  Pa,  90  percent  of  the  discharge  current  originates  from 
this  region. 

The  phenomenon  described  above  suggests  that  one  can  postulate  the  existence 
of  a  plasma  production  region  where  the  process  of  electron  emission  and 
plasma  production  are  concentrated.  To  simplify  computation,  one  can  also 
assume  that  the  plasma  parameters  are  uniform  in  this  region  and  the  position 
of  the  region  is  downstream  of  the  emitter.  Considering  the  elastic  reflec¬ 
tion  of  the  initial  electrons,  the  length  of  the  plasma  production  region  is 
chosen  to  be  Lg  =  2Xpj..  To  ensure  a  sufficient  amount  of  energy  exchange 
between  the  initial  electrons  and  other  particles,  and  to  maintain  stable 
discharge  inside  the  cathode,  the  emitter  is  designed  to  be  a  cylindrical 
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tube  whose  inner  radius  is  equal  to  or  greater  than  Xpr*  Thus,  the  length 
of  the  plasma  production  region  and  its  diameter  are  both  equal  to  2Xpj., 
and  the  effective  area  of  the  emitter  is  equal  to  Utt  Apr* 

By  using  the  ion  discharge  model  and  the  computer  program  developed  by 
PATCRS,  Wilbur  et  al.^  have  calculated  the  density  of  neutral  particles  in 
the  plasma  production  region  and  have  determined  the  relationship  between 
the  non-elastic  collision  mean  free  path  and  the  density  of  the  neutral 
particles,  as  given  by  the  following  empirical  formula: 


Ain  — 


2.83X  10^ 
«o 


-1.5; 


10%, 


(6) 


where  Gpj-  is  the  initial  electron  energy,  which  is  approximately  equal  to  the 
plasma  potential  Vp.  As  in  the  case  of  Tes  the  value  of  Vp  in  the  plasma 
production  region  does  not  vary  a  great  deal  over  a  wide  range  of  discharge 
parameter  values;  its  mean  value  is  8.7  eV. 


The  mean  free  path  of  elastic  collisions  between  the  initial  electrons  and 
Maxwell  electrons  is:^ 


Ac 
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•  rte 
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Substituting  equations  (6),  (7)  into  (5),  one  obtains: 


6.5X  lO^ttoSpr _ 

elr  2.83X  10'^-1.5no 


(8) 


Once  the  inner  radius  of  the  emitter  is  known,  the  dimensions  of  the  emitter 
can  be  determined,  and  an  appropriate  cathode  tube  (tantalum  tube)  which 
matches  the  emitter  can  be  selected. 

III.  Estimate  of  the  Effective  Emitter  Area 

The  analysis  presented  in  the  previous  section  provides  a  preliminary  esti¬ 
mate  of  the  effective  emitter  area.  From  this  estimate  one  can  determine 
whether  the  emitter  surface  is  operating  in  the  space  charge  region  or  the 
acceleration  region,  and  whether  the  effective  emitter  area  can  satisfy  the 
requirements  of  the  8-cm  thruster  (i.e..  Ip  =  0.5A). 

Analysis  shows^  that  the  emission  current  is  primarily  produced  from  thermal- 
electron  emissions;  photo-induced  emissions,  secondary  emissions  and  field 
emissions  account  for  less  than  3  percent  of  the  total  emission  current. 

The  voltage-ampere  characteristics  curve  of  the  cathode  discharge  shows  that 
the  cathode  generally  operates  in  a  relatively  flat  region  where  the  current 
is  controlled  by  the  contact  electrode  voltage;^  this  implies  that  the 
emitter  surface  may  operate  in  the  acceleration  region.  By  using  the  plasma 
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parameters  given  in  the  previous  section,  it  is  possible  to  estimate  the 
surface  field  strength  formed  by  the  sheath  layer,  and  to  calculate  its 
effect  on  the  surface  work  function.  It  can  be  shown  that  the  emitter 
surface  indeed  operates  in  the  acceleration  region. 

The  emitter  surface  field  strength  formed  by  the  plasma  sheath  layer  is: 
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where  Vp  is  the  plasma  potential,  chosen  to  be  8.7  eV; 

Aj)  is  the  Debye  length,  which  is  of  the  same  order  of  magnitude  as  the 

sheath  thickness; 

e  is  the  electron  charge; 

eo  is  the  vacuum  dielectric  constant. 

Because  of  the  presence  of  the  field  strength,  the  value  of  the  surface  work 
function  is  reduced: 
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where  is  the  effective  work  function; 

(j)g  is  the  surface  work  function  of  the  material. 

Also, 

<fi.,-<Po+aTs  (11) 


where  (f)o  is  the  work  function  of  the  material  at  a  temperature  of  absolute 
zero;  it  is  chosen  to  be  1.67  eV; 

a  is  the  temperature  coefficient  of  the  emitter  work  function;  it  is 
chosen  to  be  2.82  x  10“^. 


If  Tg  is  chosen  to  be  1273K,  then  from  equations  (11)  and  (10) ,  one  obtains 
03  =  2.03  eV  and  (|)g  =  1.83  eV.  Because  of  the  presence  of  the  field  strength 
(which  compresses  the  potential  barrier) ,  the  work  function  is  reduced  by 
0.2  eV.  Therefore,  in  calculating  the  emitter  current  density,  the  effect 
of  the  acceleration  field  on  the  work  function  should  be  taken  into  account: 


/e  =AT]  exp 


(12) 


where  je  is  the  emitter  current  density; 

A  is  the  emission  constant. 

The  discharge  current  passing  through  the  top  orifice  consists  of  two  parts: 
one  is  the  electron  current  emitted  by  the  emitter;  the  other  is  the  ion 
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current  produced  by  the  ionization  process,  i.e..  Ip  =  le  +  li-  The  current 
can  be  calculated  from  the  energy  equilibrium  relationship  in  the  plasma 
production  region,  but  is  difficult  to  accurately  compute  the  partial 
energies  of  the  different  processes  taking  place  in  the  region.  By 
incorporating  the  test  results  of  the  orificed  hollow  cathode  in  the 
calculations,  it  was  found  that  when  the  pressure  inside  the  cathode  exceeds 
400  Pa,  the  ratio  Ii/Id  is  approximately  a  constant  (~0.3)  over  a  wide  range 
of  Ip  values.  For  the  discharge-chamber  cathode.  Ip  =  0.5A,  hence 
le  =  0.35A. 

Substituting  the  values  A  =  1.2  x  10^  A/m^K^,  Tg  =  1273K,  and  (jig  =  1.83  eV 
into  equation  (12)  gives  jg  =  1.14  x  105  A/m^.  If  Xpr  is  chosen  to  be  0.5  mm, 
then  the  effective  emitter  area  of  the  cathode  is  Sr  =  27rR  x  Le  =  2TrXpr  x 
2Xpr  =  IT  X  10“6  m2,  Ig  =  Sjj  x  jg  =  0.36A,  which  meet  the  design  requirements. 

Based  on  manufacturing  and  assembly  considerations,  the  following  cathode 
dimensions  are  selected:  emitter  inner  diameter  1-1.4  mm,  outer  diameter 
2.6  mm,  tube  length  10  mm;  it  is  placed  in  a  0.2-mm-thick  tantalum  tube  with 
an  outer  diameter  of  3  mm. 

IV.  Design  of  the  Cathode  Heater^ 

The  orificed  hollow  cathode  is  a  self-sustaining  arc-discharge  thermal 
cathode.  Upon  ignition,  the  power  supply  to  the  cathode  heater  is  immediately 
cut  off  to  avoid  over-heating.  As  the  performance  of  the  cathode  emitter 
degrades  toward  the  end  of  life,  the  current  requirement  for  the  emitter  can 
no  longer  be  met  by  the  self-heating  mode;  the  energy  required  to  raise  the 
cathode  temperature  must  be  provided  by  the  heater.  At  ignition,  the  heater 
is  designed  to  heat  the  emitter  to  a  maximum  temperature  of  1200°C.  In  the 
heater  design,  precaution  must  be  taken  so  that  the  heating  filament  will  not 
evaporate  due  to  excessive  temperature,  and  the  AI2O3  layer  will  not  react 
with  the  tungsten  coating  or  the  heating  filament  to  cause  damage  to  the 
insulator  or  the  filament. 

1.  Estimate  of  the  Power  Consumption  of  the  Cathode  Heater 

As  shown  in  Figure  2,  the  cathode  heater  reaches  thermal  equilibrium  through 
heat  dissipation  in  the  form  of  radiation  and  conduction.  The  thermal 
equilibrium  equation  is: 

P//=P,+P2+P3+P4+P5+P6+P7  (13) 


where  Pjj  is  the  heating  power  of  the  heater; 

P]^  is  the  radiation  loss  through  the  top  of  the  cathode; 

P2  is  the  radiation  loss  in  the  radial  direction  through  section  A  of 
Figure  2; 

P3  is  the  radiation  loss  through  the  annular  surface  downstream  of  the 
AI2O3  layer; 

P4  is  the  radiation  loss  through  the  annular  surface  upstream  of  the 
AI2O3  layer;  P4  »  (1/2)P3; 

P5  is  the  conduction  loss  through  section  B  of  the  tantalum  tube; 
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Pg  Is  the  conduction  loss  through  section  C  of  the  tantalum  tube; 

Py  is  the  amount  of  heat  carried  away  by  the  mercury  vapor  as  it  enters 
the  cathode  and  is  heated  from  SOO^C  to  1200°C. 


Figure  2.  Temperature  Distribution  in  the  Cathode 

1.  cathode  top,  1200°C;  2.  AI2O3  layer,  1250°C;  3.  heater,  1300°C; 

4.  tantalum  tube  and  tungsten  coating,  1200°C;  5.  tantalum  tube, 

SSO^C;  6.  interface  with  mercury  vaporizer,  300°C 

As  shown  in  Figure  1,  the  heater  is  buried  in  a  layer  of  AI2O3  outside  the 
tantalum  tube;  the  emitter  is  heated  to  a  temperature  of  1200°C,  and  the 
temperature  of  the  heating  filament  is  approximately  1300®C.  The  cathode  has 
two  layers  of  heat  screens;  there  is  also  a  third  screen  which  consists  of  the 
outer  ceramic  tube  and  the  thin-wall  tantalum  sleeve.  The  temperatures  shown 
in  Figure  2  correspond  to  the  condition  when  the  emitter  is  heated  to  1200°C. 
Except  for  the  specified  temperature  value  of  1200°C  at  the  interface 
between  section  A  and  section  B  of  the  tantalum  tube,  all  other  temperatures 
are  measured  values.  In  the  calculation,  the  temperature  inside  the  dis¬ 
charge  chamber  of  the  thruster  and  the  temperatures  of  those  parts  which  are 
close  to  but  not  in  contact  with  the  cathode  are  chosen  to  be  80°C.  The 
radiated  power  is  approximated  by; 


P,= 


eaSr(T  —  To) 
n  +  1 


(14) 


where  Pj-  is  the  radiated  power; 

e  is  the  gray-body  radiation  coefficient  which  varies  with  temperature; 
a  is  the  Stefan-Boltzman  constant; 

Sr  is  the  radiation  area; 

T  is  the  temperature  of  the  radiator; 
n  is  the  number  of  heat  screens. 

The  heat  conduction  formula  is  given  by: 


P  =  liSlCP;  T]) 

L, 


(15) 


where  P^  is  the  conduction  loss; 

X  is  the  conduction  coefficient  of  the  material,  which  varies  with 
temperature; 
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St  is  the  conduction  cross-sectional  area  of  the  material; 

T2  is  the  high  temperature  of  the  conductive  material; 

Ti  is  the  low  temperature  of  the  conductive  material; 
is  the  length  of  the  conductor. 

In  equations  (14)  and  (15) ,  it  is  assumed  that  the  radiation  surface  has 
uniform  temperature  and  the  temperature  in  the  conductor  is  linearly 
distributed.  It  is  also  assumed  that  the  radiation  is  directed  only  at  the 
vertical  surface  so  that  radiation  and  conduction  losses  can  be  calculated 
independently . 

The  above  procedure  has  been  used  to  calculate  the  values  of  Pi  -  P6  for 
emitter  temperatures  of  1200®C,  1100“C,  1000®C  and  900®C;  the  results  are 
presented  in  Table  1.  The  table  also  shows  the  measured  heater  power 
corresponding  to  the  above  temperatures. 


Table  1.  Estimated  and  Measured  Heater  Power  of  the  Cathode 


\  Power  loss  through 
\  different  parts  (W) 

Calculated  value  \ 

Pi 

P2 

P3 

_ 

P4 

P5 

P6 

Ph 

(calculated) 

P 

(measured) 

Temperature  (®C) 

1200 

D 

6.17 

m 

0.8 

4.0 

16.5 

17 

1100 

0.7 

5.2 

0 

0.6 

3.8 

D 

13.4 

14 

1000 

0.5 

4.0 

0.9 

0.5 

3.8 

0.9 

11.0 

11.4 

900 

0.7 

m 

8.9 

9.2 

Py  is  neglected  because  it  is  less  than  O.IW.  The  calculated  value  of  Ph 
includes  a  multiplier  of  1.05  which  accounts  for  hot-component /cold-end 
corrections.  The  table  shows  that  the  calculated  values  are  basically  in 
agreement  with  measured  results.  Among  the  heat-loss  terms,  P2  and  P5  are 
the  major  contributors.  Also,  the  heat  filament  and  the  insulator  can 
operate  safely  even  if  the  temperature  of  the  emitter  rises  to  1300°C. 


2.  Heater  Structure 


The  heater  is  made  of 
10-turn  single-spiral 
lead -wire  on  one  end. 
0.675  Q-cm  at  1300°C. 


0.25-mm  tungsten-rhenium  wires  which  are  wound  into  a 
coil  with  an  inner  diameter  of  3.8  mm;  it  has  a  15-mm 
The  tungsten-rhenium  wire  has  a  resistivity  of 
The  length  of  the  heater  is  142  mm. 


The  resistance  of  the  heater  is  given  by; 


R  = 


dk 


(16) 
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where  R  is  the  heater  resistance  at  1300°C; 

Rl  is  the  resistivity  of  tungsten-rhenium  wire  at  1300°C; 

Lr  is  the  length  of  the  heater; 
dR  is  the  heater  diameter. 

When  the  emitter  is  heated  to  1200°C,  the  heater  current  is: 

I  =  (17) 


where  I  is  the  heater  current; 

Pr  is  the  heater  power. 

The  heater  uses  an  a.c.  power  supply  to  avoid  damage  to  the  AI2O3  layer  by 
electrolysis. 

V.  Design  of  the  Cathode  Structure 

The  cathode  components  operate  under  very  high  temperatures;  the  operating 
condition  is  dictated  by  the  operating  mode  of  the  thruster.  The  design  of 
the  cathode  structure  must  take  into  consideration  such  factors  as  the 
material  property  under  high  temperature  (including  mercury) ,  the  structural 
configuration  under  thermal  shock,  the  manufacturing  techniques,  the  thermal 
efficiency,  as  well  as  resistance  against  vibration  and  ion  bombardment  for 
the  emitter. 

The  emitter  inside  the  cathode  is  a  tubular  electrode  made  of  porous 
tungsten  soaked  into  aluminate;  it  is  produced  using  mature  technologies  and 
has  good  emission  capability  and  high  resistance  against  ion  bombardment. 

One  end  of  the  emitter  is  welded  to  the  tantalum  filament,  whose  other  end 
is  welded  to  the  inner  wall  of  the  tantalum  tube. 

As  shown  in  Figure  1,  the  outside  of  the  tantalum  tube  is  coated  with  a  layer 
of  melted  tungsten  powder  using  plasma  coating  technique  to  prevent  chemical 
reaction  between  the  AI2O3  layer  and  the  tantalum  tube  under  high  temperature. 
The  AI2O3  is  also  plasma-coated  to  provide  insulation  and  structural 
rigidity  for  the  heater.  The  coating  layer  is  rigidly  attached  to  the  base 
material.  The  coating  layer  has  a  relatively  loose  texture,  which  is 
designed  to  withstand  thermal  shock.  In  order  to  improve  thermal  efficiency, 
reduce  cathode  pre-heat  time,  and  to  avoid  cracks  in  the  coating  layer  due 
to  different  thermal  expansion  coefficients  of  the  materials,  the  parts  and 
coating  layers  all  have  very  thin  walls  and  thin  layers.  The  tungsten-powder 
layer  is  0.1-0.15  mm  thick;  the  AI2O3  insulation  layer  is  also  0.1-0.15  mm 
thick  and  its  resistance  is  several  megohms.  The  AI2O3  layer  used  to  hold 
the  heater  is  0. 2-0.3  mm  thick.  The  thicknesses  of  the  tantalum-foil  heat 
screen,  the  tantalum  tube,  and  the  ceramic  tube  are  respectively  0.01  mm, 

0.2  mm,  and  0.4  mm.  The  thln-plate  cathode  top  with  a  center  orifice  is 
made  of  porous  tungsten  material;  this  material  is  chosen  over  dense  tungsten, 
molybdenum  or  tantalum  in  order  to  keep  the  orifice  from  expansion,  distortion 
or  fusion.  It  is  also  easily  welded  to  the  tantalum  tube  using  electron-beam 
welding  techniques;  the  seam  is  smooth,  free  of  cracks,  and  will  not 
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deteriorate  under  ion  bombardment.  The  heat  filament  is  made  of  tungsten- 
rhenium  material  which  has  a  higher  re-crystallization  temperature  and  better 
ductility  and  malleability  than  pure  tungsten;  it  also  has  good  resistance 
against  thermal  shock  and  vibration. 

Thermal-environment  simulation  tests  show  that  in  the  temperature  range  80°C- 
1030°C,  the  cathode  can  withstand  13,000  thermal  cycles  each  consisting  of 
4-min  temperature  rise  and  16-min  temperature  drop.  In  the  temperature  range 
1050“C-1100®C  with  no  discharge,  the  cathode  can  operate  for  30,000  hours 
without  damage,  and  all  the  cathode  and  heater  parameters  remain  stable.  The 
heater  resistance  shows  no  apparent  changes  and  the  resistance  of  the  AI2O3 
insulation  layer  remains  at  several  thousand  ohms. 
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[Text]  Abstract 

The  self-focusing  process  and  self-focusing  ability  of  a  double  helical 
wlggler,  commonly  used  in  a  Raman  free  electron  laser  (FEL) ,  is  analyzed. 

The  conditions  for  steady  beam  transportation  are  obtained  and  compared  to 
the  experimental  beam  data  from  a  Raman  FEL  without  axial  magnetic  field.  In 
addition,  the  FEL  radiation  data  obtained  on  the  EPA-74  electron  pulse  accel¬ 
erator  is  analyzed.  The  theoretical  analysis  is  in  good  agreement  with  the 
experimental  results. 

PACC:  4255T,  4060,  7845. 

I.  Introduction 

In  an  FEL,  the  steady  transport  of  the  electron  beam  is  influenced  by  the 
wiggler  field  to  produce  transverse  oscillation.  The  coupling  between  the 
radiation  field  and  wlggler  field  creates  a  ponderomotive  potential  which 
captures  electrons.  This  leads  the  bunching  of  electrons.  The  kinetic 
energy  lost  by  the  bunched  electron  beam  is  converted  into  electromagnetic 
radiation  which  coherently  amplifies  the  radiation  field.  Consequently,  a 
high-power,  high-efficiency  coherent  radiation  output  is  produced. 

FEL's  may  be  divided  into  Raman  and  Compton  FEL's.  Regardless  of  the  type, 
in  order  to  understand  the  physics  Involved,  it  is  necessary  to  investigate 
the  state  of  electron  motion  and  the  beam  transportation  process.  Freund^ 
and  Fajans^  extended  the  work  done  by  Friedland^  to  obtain  some  semi- 
empirical  results  on  the  state  of  electron  motion  in  the  magnetic  field  of  a 
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double  helical  wiggler  without  taking  space  charge  into  account.  These 
results  have  an  important  impact  on  the  understanding  of  the  physics 
associated  with  an  FEL.  Furthermore,  both  theoretical  and  experimental 
studies  show  that  the  presence  of  an  external  guide  field  could  enhance  the 
power  and  efficiency  of  an  FEL  under  certain  conditions.  Nevertheless,  this 
raises  issues  such  as  difficulty  of  identifying  the  radiation  mechanism.^ 

In  order  to  overcome  these  difficulties,  we  conducted  a  theoretical  and 
experimental  study  on  the  Raman  FEL  without  an  axial  magnetic  field.  In 
this  experiment,  from  a  280-A  electron  beam  at  560  keV,  the  radiation  output 
obtained  has  a  peak  power  of  over  7.6  MW  and  a  frequency  of  35-38  GHz. 

II.  Analysis  of  the  Electronic  State  of  Motion 

Let  us  assume  that  the  field  experienced  by  an  electron  after  an  electron 
beam  enters  the  homogeneous  zone  of  the  wiggler  magnetic  field  through  its 
adiabatic  inlet  is 


5  = 


(1) 


where  the  subscripts  ex  and  se  represent  the  external  field  and  self-field, 
respectively.  When  analyzing  the  steady-state  electron  motion,  the  smaller 
radiation  field  is  neglected.  Furthermore,  the  higher  harmonics  of  the 
wiggler  field  are  omitted.  Then,  the  external  magnetic  field  may  be 
expressed  as^ 

Efx  ^ I co&Xffif^  ~  id  */ sin  sin XqC, ]  (2) 

where  Xq  =  6o  ~  Xq  =  k(jjro.  The  static  self-field  of  the  beam  is 


ii.e —  (l/2e|,)»ofroe,„.  (3) 

Here,  it  is  assumed  that  the  beam  current  has  axial  symmetry  and  its  density 
is  no . 

Let  us  consider  steady-state  electron  motion  and  assume  that  its  kinetic 
energy  remains  unchanged  during  transportation.  Then,  the  equation  of 
motion  is 


dvjdt  —  —eCE  +  Co  X  B')jrmo.  (4) 

Under  the  influence  of  the  field  given  in  equations  (2)  and  (3) ,  the  steady- 
state  motion  of  the  electrons  near  the  axis  (Xo  <<  1,  i.e.,  rg  s  0)  is  helical 
with  a  period  that  corresponds  to  the  period  of  the  wiggler. 5  The  focusing 
force  is  from  the  presence  of  an  axial  gradient  of  the  wiggler  field.  This 
effect  may  be  recognized  from  the  e^o  component  in  equation  (4) ; 
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ro/(2rrJ), 


(5) 


d^rjdt^^"  —  +  t/,  sinX,)/rmo  +  ®p  * 

where  col  —  is  the  non-relativistic  plasma  frequency.  Equation  (5) 

also  reveals  that  in  order  to  have  a  stable  transportation  of  the  beam 
current,  we  must  satisfy  <sin  Xo>  >  0.  The  symbol  <>  represents  the  average 
over  a  wiggler  period.  Hence,  only  Xo  is  a  constant.  This  means  that  the 
helical  motion  of  the  electron  should  be  synchronous  with  the  wiggler  field. 
For  simplicity,  we  substitute  the  velocity  without  taking  space-charge 

effect  into  account,  into  equation  (5).^  Therefore,  the  conditions  for 
stable  transportation  beam  current  are: 

<  0.369,  (6) 

Qi>colK^rrl),  (7) 


where  .0,  ■—  ^ 


However,  electrons  far  away  from  the  axis  (not  meeting  the  condition  Xo  1). 
such  as  the  electrons  on  the  edges  of  the  beam  used  in  reference  6  and  in 
this  experiment,  cannot  use  the  above  approximation.  In  this  case,  a  guide 
center  approximation  methodl»2  may  be  used  to  perform  the  analysis. 


As  shown  in  Figure  1,  electron  motion  is  decomposed  into  two  parts;  the  slow 
motion  at  the  guide  center  and  the  fast  motion  around  the  center;  i.e.. 


Vo  =  V  +  Vg,  Xo  = 


Xg  +  r. 


Moreover,  |vg|  «  |v|  and  |rg|  »  [r 


Here,  the 

subscript  g  represents  quantities  relevant  to  the  guide“center .  This  physical 
picture  is  basically  consistent  with  numerical  simulation.^  Hence,  the  field 
exerted  upon  an  electron  can  be  expressed  as 


fi(r,)  =  iB(r,)  +  (r  •  V,)B(r,), 

£:(ro)  =  +  (r  •  V,)£(r,),  (8) 


where  V,  =  d Idr,,  By  substituting  (8)  into  (4)  and  taking  the  average  over  a 
wiggler  period,  we  have 


dVJdt  --i-  [£(rj  +  {V,  X  (r  •  V,)£(r,)> 

rwj 

+  <K  X  (£(r,)  +  (r  •  V,)£(r,))>] 
dVIdt  +  +  (r  •  V,)B(r.) 

-  <(r  .  V,)B(r,)»  +  K  X  (£(0  +  (r  .  V,)£(0) 

-<KX  (jB(r,)  +  (r.  V,)£(r,))>],  (10) 

Here  the  S3nnbol  <>  represents  the  average  over  a  wiggler  field  period.  The 
focusing  process  of  the  wiggler  field  may  be  analyzed  and  the  state  of 
electron  motion  solved  based  on  equations  (9)  and  (10).  First,  let  us  solve 
the  fast  component  of  electron  motion.  Equation  (10)  is  too  complicated  to 
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solve  accurately.  Accurate  to  the  first  order  of  l(=k(jjr),  equation  (10)  may 
be  written  as 

dvjdt  ^  u)lrl(2rr\')e,—  eV  X  iBir,)  +  (r  •  Vi)S(r,))/(r;wo).  (11) 

In  order  to  analyze  the  steady-state  motion  of  the  beam  and  its  focusing 
process,  it  is  necessary  to  analyze  the  component  of  equation  (11); 

f  —  (.Vaytr  ■=  Mp  r/(2rri  )«r  —  •  [79/1(1^)  sin  X, 

+  F, 7,(1^)  sin  X/ij  4- F,/2(i,)cosX,  sin(0  — 

+  (T^a/i(l<)?'/>'*)sinX  +  —  d^)sinX, 

+  (V.^2(Zr)/r,)  sin  (X  +  0  -  &,).  (12) 

According  to  our  assumption,  <sin  Xg>  ~  0,  <sin(6-0g)>  »  0,  <cos(6-0g)>  »  0, 
and  <sin(X+0-0g) >  «  0  in  the  electron  motion.  Hence,  these  corresponding 
terms  do  not  have  any  focusing  effect.  Only  when  <sinX>  >  0  does  the 
gradient  of  the  axial  component  of  the  wiggler  field  have  a  constraining 
effect  on  the  fast  motion  of  the  electrons.  In  this  case,  X  is  a  constant. 
In  order  to  find  a  specific  solution  of  this  fast  motion  component,  as  an 
approximation  (as  well  as  for  simplicity) ,  the  space-charge  effect  of  the 
beam  current  and  the  gradient  component  of  the  field  are  not  considered. 
Hence, 


dVjdt  Rtf  —e(y  X  (13) 

Let  us  choose  the  following  helical  coordinates;  eq  =  erCosX  -  e0sinX, 

®2  =  CrslnX  +  e0CosX,  €3  =  Then,  the  steady-state  motion  of  the  electrons 

with  a  period  k^^z  can  be  solved  as  follows: 

V- F,e.  +  F„ei,  (14) 

where  V^j  =  (Xe) /k;^,  =  eB^j/ymo,  and  X  =  ir/Z.  Therefore,  the  major  fast 

motion  component  is  a  helical  motion  around  the  guide  center. 


Figure  1.  Beam-Current  Cross  Section 
and  Selection  of  Coordinate  System 
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Nevertheless,  equation  (14)  does  not  take  into  account  the  effect  of  limited 
helical  radius  on  the  motion.  When  the  wiggler  field  is  relatively  strong, 
this  effect  cannot  be  neglected.  A  semi-empirical  expression  which  takes 
into  consideration  the  limited  helical  radius  is  introduced  in  the  analysis. 

-  2QJ,iX,^hix'))).  (15) 

Substituting  Vq  into  (12) ,  one  gets  the  stability  conditions  for  fast  elec¬ 
tron  motion  as  follows: 


(.QJMY  >  (ol  /(4rrO,  (16) 

<  0.369.  (17) 

Substituting  (15)  into  (9),  it  is  possible  to  obtain  the  equation  of  motion 
for  the  guide  center. 


dV dt  oj^e,  X  Vg  +  Ato/2r,  —  V,^col^crilYcej,  (18) 

where  w, - Q^V^jV,,  - Q^k.u.V^  +  mll2rr\ , 

The  transverse  component  may  be  written  as 


U^ldi-  +  Qf)id^ldt^  +  (  y  )  -  0 .  (19) 

where  Q\,i  =  (wj  —  Aco)/2±m,(w^  —  2Aw)‘^V2. 

When  the  wiggler  field  amplitude  is  small,  we  get 

Q\.i  =  (-Qi  —  /rri)/2.  (20) 

Hence,  equation  (19)  becomes 

(^V*^  +  PO(^*)  =  0.  (21) 

The  solution  to  equation  (21)  is  the  transverse  motion  of  the  electron  guide 
center,  i.e..  Betatron  motion.  In  order  to  have  stable  electron  motion,  in 
addition  to  having  a  steady  fast  helical  motion,  its  Betatron  motion  must 
also  be  stable,  i.e.. 


Qi>(ollrr\.  (22) 

The  beam- exp ending  effect  from  the  space  charge  is  offset  by  the  presence  of 
an  axial  component  of  the  wiggler  field  (when  Xg  ^  1,  this  component  is  pro¬ 
portional  to  rg) . 
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In  summary,  the  transverse  inhomogeneity  of  the  double  helical  wiggler  can 
confine  the  electron  beam.  However,  the  beam-confining  process  is  far  more 
complicated  than  that  described  in  reference  7.  In  order  to  achieve  steady 
beam  current  transportation,  the  relevant  parameters  must  satisfy  equations 
(17)  and  (22). 

In  addition,  it  is  revealed  in  reference  2  that  when  space  charge  is  not 
considered  and  if  the  axial  magnetic  field  is  small  or  absent,  there  is  a 
substantial  difference  between  the  numerical  simulation  and  the  approximation 
results  at  the  center  due  to  faster  Betatron  motion.  However,  the  analysis 
shown  above  shows  that  Betatron  motion  can  be  slowed  down  by  space-charge 
effect  and  the  Raman  FEL  theory  can  still  be  applied  when  there  is  no  axial 
magnetic  field. 

III.  Theoretical  Considerations  of  Raman  FEL  Without  Axial  Magnetic  Field 

The  basic  physical  process  in  a  Raman  FEL  can  be  viewed  as  an  energy  transfer 
process  involving  the  coupling  of  the  wavegulded  electromagnetic  wave  and  the 
space  charge  wave  to  amplify  the  electromagnetic  wave^  in  an  applied  wiggler 
field.  The  dispersion  relations  are: 

£0^  ■=  c'^lC  +  w?,  (23) 

w  =  (^  +  +  pico^/r,r^'\  (24) 

where  +  pital  Jr  is  the  effective  cutoff  frequency  of  the  waveguide, 

Wco  is  the  cutoff  frequency  of  the  waveguide  in  vacuum  and  p2  can  be  treated 
as  a  filling  factor  for  the  energy  of  the  electromagnetic  wave  in  approxima¬ 
tion.  The  energy  transfer  primarily  takes  place  at  the  frequency  synchronous 
to  the  slow  space-charge  wave,  i.e., 

(25) 


where  —  piWp/(r,r"^^/S,c)  is  the  effective  wave  number  of  the  wiggler 

field.  Theoretical  analysis^  and  experimental  results^O  indicate  that  the 
low-frequency  radiation  gain  is  low  for  equation  (25) .  Furthermore,  it  is 
not  of  concern  to  this  work,  and  therefore  will  not  be  taken  into  considera¬ 
tion. 

The  Raman  FEL  gain  is  the  ratio  of  the  radiation  power  at  z  to  that  at  1 
(<z)ll  (relative  to  the  amplifier  1=0); 


G=  <.(^)/K0=  loC>r)/a(/)|%  (26) 

where  la(z)/a(l)|  is  the  ratio  of  the  electric  field  amplitudes,  which  can  be 
expressed  as 


o(s)/a(/)  ==•  5]  —  0)> 


(27) 


27 


The  summation  is  with  respect  to  all  mutually  interactive  radiation  fields, 
Ai  is  the  normalized  amplitude  of  a  radiation  field,  and  is  its  complex 
propagation  constant.  They  can  be  derived  from  the  poles  and  remainders  in 

dispersion  equations 


A(sO  =  D/F, 


(28) 


where  D 6  +  +  id.Xsk- 6  -  6^  +  W,),  F  -  6  +  6,  -h  id,)  ^ 

•  —  6  —  6^  +  id,)  +  Q,  5  ■=  is  an  "out-of-tune  factor"  of  the  FEL, 

Q  ““  PifixkiTidlil  +  fi,yi(A^p\)  represents  the  intensity  of  energy  exchange  between 
the  two  waves,  is  the  normalized  plasma  frequency, 

toAr,exp(— 2^VA)(fi/A)''/(ri/Ji^)  is  the  normalized  Landau  damping  rate. 
Here,  Q  pico^rJicar^^),  A  =  Ar,/(ri/?i).  Therefore,  for  a  "cold"  electron  beam, 
Ar,/r, ->0,  d,->0.  Or,  hy^Hz  “  3%  is  a  constant;  when  the  beam  current 
increases  to  a  certain  value,  0^  ->  0.  This  indicates  that  Landau  damping  is 
very  small  and  can  be  neglected.  In  our  experimental  system,  Ayz/Yz  “  3%. 

As  long  as  I-p  >.  40A,  the  effect  of  Landau  damping  can  be  neglected  and  the 
beam  can  be  considered  as  a  cold  beam.  Thus,  equation  (28)  can  be  written  as: 

A(i80  =  (8^-6  +  e,,)(8k  —  d  —  d^)fC8kCs^  -d  +  d^)Csi^  —  d 

+  (29) 

When  0  =  -0p  (which  corresponds  to  the  situation  of  (25)),  the  maximum  gain 
rate  (t  =  Imk)  of  the  Raman  FEL  is 

I  r|  -=  yy  {p^r.r^Mn.  (30) 

This  is  the  space  gain  rate  in  the  linear  gain  region.  When  most  electrons 
are  captured  by  the  ponderomotive  force  and  transfer  some  energy  to  the 
radiation  field,  the  radiation  output  saturates.  This  corresponds  to  the 
situation  in  our  experiment  when  the  mutual  interaction  zone  was  approxi¬ 
mately  1  m  long.  The  maximum  radiation  output  efficiency  is  determined  by 
that  state.  However,  in  principle,  linear  theory  should  not  be  applied  to 
theoretically  analyze  FEL  efficiency.  Instead,  a  nonlinear  theory  that 
describes  the  saturated  state  should  be  employed.  Nevertheless,  theoretical 
analysis  shows  that  linear  theory  may  be  used  to  estimate  FEL  eff iclency . 
Furthermore,  for  an  FEL  with  a  magnetic  wlggler,  the  estimated  value  based  on 
linear  theory  is  in  good  agreement  with  the  result  from  self-consistent  non¬ 
linear  theory. The  saturation  efficiency  of  FEL  isrH 

v-(:fiJ^)''^rXl+r)ip,a>^Kcnrn.  (31) 

IV.  Comparison  of  Experimental  Results  to  Theoretical  Analysis 

The  Raman  FEL  experimental  system  without  an  axial  magnetic  field  is  shown  in 
Figure  2.  This  system  was  successfully  used  to  carry  out  spontaneous  radia¬ 
tion  amplifier  Raman  FEL  experiments. 
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Figure  2.  The  EPA-74  Raman  FEL  System  Without  an  Axial  Magnetic 
Field 

1.  cathode;  2.  lens;  3.  wiggler;  4.  drift  tube;  5.  output 
horn;  6.  attenuator;  7.  dispersion  line;  8.  band-pass 
filter;  9.  crystal  rectifier;  10.  oscilloscope 

An  electron  beam,  560  keV  in  energy,  100-300  A  in  intensity,  3  percent  in 
energy  dispersion,  0.06  cm-rad(ei^)  in  emissivity,  and  approximately  1  cm  in 
diameter,  enters  the  1.7-m-long  drift  tube  which  is  pumped  to  vacuum.  The 
diameter  of  the  tube  is  1.6  cm.  A  circularly  polarized  wiggler  field,  3.45  cm 
per  period,  is  generated  by  a  1.5-m-long  concentric  double  helix  mounted  on  a 
nylon  rack.  A  7-cycle  adiabatic  input  segment  guarantees  the  adiabatic 
transition  of  the  beam  to  the  homogeneous  zone  of  the  wiggler  field.  It  has 
a  helical  trajectory  such  as  that  described  in  section  II  of  this  paper. 

When  the  wiggler  field  was  measured  to  be  1.25  kGs  by  using  an  8-mm  disper¬ 
sion  line  measuring  system^^  and  a  microwave  calorimeter,  the  radiation 
frequency  is  35-38  GHz,  energy  is  41-48  mJ,  and  peak  power  is  7.6  MW.  The 
following  is  an  analysis  of  the  experimental  data  based  on  the  results 
described  in  sections  II  and  III. 

1.  Beam  Current  Transportation 

Figure  3  shows  the  beam  current  transportation  results  without  an  axial 
magnetic  field.  It  shows  that,  within  a  certain  range,  the  transported  beam 
current  increases  with  rising  wiggler  field.  However,  when  the  wiggler  field 
amplitude  is  too  large,  beam  current  drops.  A  similar  result  was  obtained  at 
the  Naval  Research  Laboratory, ^5  This  result  was  analyzed  using  equations 
(17)  and  (22) .  From  (17) ,  in  order  to  maintain  the  desired  helical  motion, 
the  wiggler  field  must  not  be  too  high.  For  our  experimental  system,  it 
should  be  <  1700  Gs.  When  the  wiggler  field  reaches  a  certain  value  (or 
the  beam  diameter  is  enlarged  due  to  emissivity  or  space-charge  effect  and 
the  magnetic  field  at  the  edge  of  the  beam  is  increased  to  a  certain  value) , 
equation  (17)  is  no  longer  satisfied  and  the  electron  beam  suffers  some  loss. 
Hence,  in  Figure  3(a)  when  the  wiggler  field  reaches  a  certain  level,  the 
beam  current  declines  with  Increasing  magnetic  field.  With  increasing  trans¬ 
portation  distance,  the  beam  current  also  falls  due  to  losses,  as  shown  in 
Figure  3(b).  Curve  3  in  Figure  3(b)  was  obtained  after  the  coaxial  system 
was  fine-tuned.  Therefore,  the  induced  current  and  transported  current  are 
larger  than  those  in  curves  1  and  2. 
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(a)  transported  current  vs, 
wiggler  field  amplitude 
(at  z  =  70  cm) 


0  20  40  60  80  100  120 


z(cm) 

(b)  transported  current  vs. 
transmission  range 


Figure  3 


In  order  to  have  stable  transportation  of  beam  current,  the  Betatron  motion 
must  also  be  stable.  Equations  (17)  and  (22)  can  be  used  to  obtain  the 
threshold  stable  beam  current  as  shown  in  Figure  4.  One  can  see  that  the 
stronger  the  wiggler  field  is,  the  higher  is  the  current  that  is  allowed  to 
maintain  a  stable  beam  motion.  This  is  consistent  with  the  pattern  observed 
in  the  experiments.  The  numerical  difference  is  caused  by  these  facts; 

1)  a  large  number  of  electrons  are  lost  in  the  adiabatic  transition  segment 
of  the  wiggler  field  (electrons  that  passed  the  segment  should  obey  the 
pattern  shown  in  Figure  4  even  in  the  absence  of  other  fields)  and  the  beam 
current  at  the  start  of  the  homogeneous  zone  of  the  wiggler  field  is  smaller; 

2)  the  injected  beam  is  not  an  ideal  parallel  beam,  which  leads  to  divergence 
and  losses;  and  3)  equation  (22)  is  derived  under  certain  ideal  conditions. 
This  is  somewhat  different  from  reality,  especially  when  the  wiggler  field 
amplitude  is  too  strong.  However,  the  above  analysis  shows  that  as  long  as 
we  can  deliver  a  high-quality  electron  beam  below  the  threshold  to  the 
homogeneous  zone  of  the  wiggler  field,  the  wiggler  field  is  capable  of  con¬ 
fining  the  beam.  Moreover,  this  confining  force  is  very  strong.  Kirkapatric^ 
used  a  double  helix  wiggler  field  with  a  3.14-cm  period  to  confine  and  trans¬ 
port  a  2.2-MeV,  930-A  electron  beam  across  the  wiggler  field  of  1.3-1. 5  kGs. 


Figure  4.  Beam  Current  Confining  Threshold, 

Y  =  2.067,  =  3.45  cm 

Curves  1,  2  and  3  correspond  to  experiments  1, 

2  and  3.  Distance  z  starts  from  the  entrance 
of  the  wiggler.  B^j^  =  1.25  kC  for  curves  1  and 

3  and  B^j  =  1.44  kC  for  curve  2. 
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2.  Analysis  of  Radiation  Output 

In  these  experiments,  because  the  beam  radius  is  relatively  large,  the  three- 
dimensional  theory  described  in  sections  II  and  III  is  used.  Equations  (15) 
and  (25)  can  be  used  to  determine  the  radiation  frequency  and  mode  of  the 
Raman  FEL  as  a  function  of  other  parameters,  as  shown  in  Figure  5.  At  = 

3.45  cm  and  a  beam  energy  of  560  keV,  in  our  experimental  system  only  the 
TE]^2.  "^^01  may  exist.  There  are  two  reasons  why  it  is  hard  to 

operate  in  mode  TMqi:  1)  In  the  experiment,  in  order  to  ensure  good  beam 
transportation  to  conduct  an  effective  Raman  FEL  experiment,  is  chosen  to 

be  between  1  and  1.5  kGs.  This  kind  of  magnetic  field  falls  outside  mode 
TMqi.  2)  It  is  difficult  for  an  electron  beam  injected  concentrically  with 
respect  to  the  drift  tube  with  current  concentrated  near  the  axis  to  excite 
mode  TMqi,  whose  transverse  electric  field  component  on  the  axis  is  zero. 
Figure  6  shows  typical  beam  current,  beam  voltage  and  radiation  output  wave¬ 
forms.  The  narrowing  of  the  radiation  output  pulse  might  be  caused  by  the 
narrowing  or  variation  of  the  beam  current  and  energy. 


(b) 


Figure  5.  Radiation  Output  Frequency  vs.  (a)  Wiggler  Field 

(b)  Electron-Beam  Energy.  I  is  the  experimental  value. 


(a)  Electron  beam, 
voltage  waveform 


(b)  Transported  (c)  Radiation  output 

current  waveform  waveform 

(z  =  70  cm) 


Figure  6.  Typical  Signal  Waveforms 
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Under  the  conditions  that  B^j  is  1.25  kGs,  beam  current  is  280  A,  beam  energy 
is  560  keV,  and  beam  diameter  is  1  cm,  the  measured  FEL  gain  curve  is  as 
shown  in  Figure  7.  Under  these  conditions,  «  0,85  and  p2  ~  0.6. 

Figure  7(b)  shows  the  spatial  radiation  growth  rate  curve  derived  from 
equation  (30) .  It  is  in  excellent  agreement  with  measured  results . 


a 

a: 


(a)  Radiated  power  vs. 
wiggler  field  length 


(b)  Growth  rate  t  vs. 

dotted  line:  theoretical; 
+:  experimental  point 


Figure  7 

After  passing  through  the  linear  gain  region,  radiation  saturates  at  approxi¬ 
mately  1  m  in  the  mutual  Interaction  zone.  The  saturated  power  is  7.6  MW, 
which  corresponds  to  a  5  percent  efficiency.  This  is  quite  different  from 
the  theoretical  10  percent  efficiency,  probably  because  of  beam  quality  and 
precision  of  the  wiggler  field. 

V.  Conclusions 

From  the  above,  we  can  conclude  that: 

(1)  Due  to  inhomogeneous  spatial  field  distribution,  electrons  away  from  the 
axis,  in  addition  to  the  fast  helical  motion  around  the  guide  center  at  a 
period  of  also  have  a  relatively  slow  motion  with  respect  to  the  guide 

center.  The  space-charge  effect  further  slows  down  this  motion.  Conse¬ 
quently,  within  a  certain  range,  the  electron  motion  without  an  axial 
magnetic  field  can  still  be  adequately  described  by  the  guide  center  approxi¬ 
mation  in  the  analysis  of  the  physical  processes  of  a  Raman  FEL. 
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(2)  Experimental  results  obtained  with  the  EPA-74  Raman  FEL  agree  well  with 
the  theoretical  analysis. 
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[Text]  Abstract 

In  this  paper  the  advantages  of  the  realisation  of  a  far  infrared  (FIR)  free 
electron  laser  (FEL)  are  shown:  it  requires  an  electron  beam  with  less 
stringent  characteristics,  less  expensive  and  it  is  unique  source  for  non¬ 
linearity  and  coherence  in  FIR.  The  dependence  of  the  FEL  performance  on 
parameters  of  the  electron  beam  and  undulator  is  described.  For  a  given 
beam  emittance  there  is  a  specific  wavelength  below  which  the  gain  rapidly 
becomes  very  low.  The  important  physical  phenomena  in  the  FIR  and  possible 
uses  of  FIR  FEL  are  given  in  broad  outline,  these  studies  could  be  important 
for  the  possible  development  of  a  THz  electronics. 

Key  words:  Far  infrared;  Free  electron  laser. 

I.  Introduction 

The  free  electron  laser  (FEL) ,  based  on  stimulated  emission  of  radiation  from 
ultrarelativistlc  electron  beams  in  an  undulator  (periodic  transverse 
magnetic  field)  has  caused  great  interest  for  its  possibilities  to  extend  the 
range  of  electron  beam  devices  (travelling  wave  tubes  etc.)  to  shorter  wave¬ 
lengths  than  the  microwaves  (cm-waves)  used  up  to  now.  Since  the  first 
operation  of  a  FEL  oscillator  in  1976,  although  slowly  (as  their  construction 
is  not  easy),  a  number  of  FEL  oscillators  have  been  operated,  in  the  range 
from  mm-waves  to  the  visible  range.  Some  of  them  are  beginning  or  will  soon 
be  used  as  user  facilities . Efforts  are  under  way  to  extend  the  range  of 
FELs  as  far  as  possible  into  the  VUV;  however  their  possible  actual  use  is 
not  in  the  near  future. 

The  far  infrared  (FIR)  region  of  the  spectrum  is  the  range  extending  approxi¬ 
mately  from  40  ym  to  1  mm  wavelength. 
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In  this  paper  we  want  to  show  the  advantages  of  the  realisation  of  a  far 
infrared  (FIR)  FEL  particularly  for  three  reasons: 

1)  The  operation  at  longer  wavelengths  is  less  delicate  and  requires  an 
electron  beam  with  less  stringent  characteristics. 

2)  As  it  requires  a  lower  energy  machine,  a  FIR  FEL  is  less  expensive  than 
in  the  near  IR. 

3)  There  is  a  lot  of  interesting  physics  to  do  in  the  FIR  with  an  intense, 
short  pulse  and  tunable  source,  and  the  FEL  is  the  only  source  having  all 
these  characteristics  in  the  FIR. 

First  we  describe  the  dependence  of  the  FEL  performance  on  its  parameters, 
and  then  we  indicate  some  of  the  fields  of  application  in  the  FIR. 

II.  FEL  Characteristics 

Let  ymc^  be  the  electron  energy  (y  »  1)  and  I  the  peak  current,  Xo  the 
undulator  period,  and  K  =  eBoXo/2'rrmc  its  "deflection  parameter".  The  range 
of  operating  wavelengths  is  defined  by5»6 

(1) 

The  efficiency  can  be  of  the  order  of  1/2  N  (N  =  number  of  period) ,  but  even 
higher  with  a  tapered  undulator.  The  peak  power  is  then  Pl  '=‘  y™c2i/2N. 

The  important  parameter  determining  how  easily  oscillation  can  be  obtained  is 
gain.  The  dependence  of  gain  on  wavelength  is  of  course  depending  on  which 
parameters  we  keep  constant. 

A  condition  that  practically  has  to  be  satisfied  to  have  a  useful  gain  is  that 
the  electron  beam  is  matched  to  the  (diffraction-limited)  mode  of  the  emitted 
radiation  in  the  cavity.  That  means  that  the  electron  beam  transverse  sizes 
are 


O'* — O'*  —  (iV2w) 


(2) 


(L  =  NXo  length  of  the  undulator)  and  the  divergences 

(3) 

In  accelerator  terms,  this  means  emittance 

s=^V2»f,  (^) 

and  full  coupling  (equal  vertical  and  horizontal  emittances) ,  and  beta  func¬ 
tion  of  the  order  of  L/2ir.  A  smaller  emittance  would  give  the  same  gain, 
but  it  would  be  more  difficult  to  get  the  same  current. 
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Another  necessary  condition  is  that  the  electron  energy  spread  should  be 
small : 


/iy/y<l/2N,  (5) 

The  gain  is  also  reduced  if  the  electron  bunch  is  shorter  than  NX;  so  we 
should  have 

At>‘N'Klc,  (6) 

At  is  the  sliping  time  of  one  optical  wavelength  for  each  wiggler  period,  then 
it  might  be  useful  to  use  a  "debuncher"  (or  "energy  compressor")  where  Ay  is 
decreased  while  At  is  increased.^ 

In  case,  changing  wavelength,  we  change  the  beam  size  and  divergence  in  order 
to  always  satisfy  the  previous  conditions,  the  gain  can  be  expressed  in  the 
form 


gc^5.9-I-Nyy,  (7) 

JLa 

where  1^,  -  17045  A  is  the  Alfvdn  current,  I  is  the  peak  current. 

More  conservatively,  taking  into  account  that  When  in  eqs.  4,  5  equalities 
hold,  there  is  a  gain  reduction  of  a  factor  ~2,  we  could  put  1  instead  of  5.9 
in  eq .  7 . 

If  the  gain  is  high,®>^»^  the  single-pass  gain  increases  exponentially: 

gH^~exp(2.7g^^^),  (8) 

and  there  is  the  large  accompanying  optical  phase  shift. 

But  eqs.  7-8  are  applicable  only  in  the  conditions  (4,  5),  while  when  some 
value  exceeds  these,  the  gain  decreases  very  rapidly.  Then  for  a  given  beam 
emittance  there  is  a  wavelength  below  which  the  gain  rapidly  becomes  very  low. 

From  the  point  of  view  of  designing  a  machine  for  a  given  output  wavelength, 
its  emittance  can  grow  in  proportion  to  X.  If  the  obtainable  current  can  be 
estimated  to  be  proportional  (for  a  given  duty  cycle  and  bunch  length)  to 

ta^B.y^e^v,  (9) 

with  Bg  depends  on  the  electron  gun  and  is  of  the  order  of  lO^-lO^^  A/m^ 
(possibly  up  to  10^2) . 

In  practice  then  eq.  7  can  be  written  as: 


g  =  AN^ 


(10) 
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where  A  ~  150  m“2  for  Be  r  10®  A/m^  if  we  use  the  factor  1  (instead  of  5.9) 
in  eq .  7 . 

We  see  that  even  in  the  range  of  wavelength  longer  than  the  minimum  defined 
by  the  beam  emittance,  the  gain  tends  to  increase  considerably  with  wave¬ 
length. 

Also,  the  lower  value  of  y  necessary  to  reach  a  longer  X  make  the  machine  less 
expensive. 10 

A  table  of  examples,  for  an  undulator  with  Xo  =  5  cm,  N  =  75,  K  =  1  and  an 
electron  gun  with  Bg  ~  10®  A/nfi: 


y 

60 

40 

20 

I  (A) 

1 

2.2 

9 

X  (pm) 

10.4 

23.4 

94 

g 

0.5% 

2% 

15% 

The  effect  of  energy  spread  does  not  depend  on  wavelength  (see  eq.  5):  but 
the  availability  of  a  higher  gain  allows  to  use  a  lower  value  of  N,  a  better 
tolerance  to  hy/y  can  be  obtained,  at  the  expense  of  a  decrease  in  gain  by  a 
factor  n2;  but  the  decrease  in  gain  due  to  energy  spread  is  exponential,  so 
an  optimum  is  found  for  a  value  comparable  with  the  condition  (5) ,  with  a 
peak  current  proportional  to  1/N.  An  "energy  compression"  could  be  good  both 
for  satisfying  eq.  5  and  eq.  6  (at  the  expense  of  a  reduction  in  peak 
current) . 

The  temporal  structure  of  the  light  pulse  (as  long  as  the  beam  length  is 
>  NX)  is  equal  to  that  of  the  electrons,  e.g.,  in  linacs  this  is  a  few 
picoseconds . H 

The  linewidth  is  close  to  the  Fourier  transform  limit: 

(11) 

III.  Other  FIR  Sources 

1)  Up  to  now,  most  of  the  work  on  FIR  has  been  done  with  "blackbody"  sources 
with  Fourier-transform  spectroscopic  detection  (with  a  Michelson 
interferometer) .  These  are  very  weak  incoherent  sources  and  their  use  is 
limited  to  absorption  spectroscopy. 

2)  Information  on  this  spectral  region  can  be  obtained  by  Raman  spectroscopy 
in  the  visible  range.  The  physical  information  obtained  is  often  comple¬ 
mentary  to  the  previous  measurements,  as  there  are  different  selection  rules. 

3)  Several  laser  sources  have  been  developed  in  the  FIR:  methanal,  formic 
acid,  etc.,  pumped  by  CO2  lasers.  These  are  coherent  continuous-wave  sources 
of  many  mW,l®  and  some  studies  of  coherent  radiation-matter  interactions  have 
been  done  with  them.^^  However,  although  there  are  many  lines,  they  are 
almost  not  tunable. 
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4)  A  widely  tunable  continuous-wave  FIR  source  can  be  obtained  by  mixing  two 
singlemode  CO2  lasers  on  a  point-contact  diode. 15  However,  the  power  is  in 
the  nano  watt  (nW)  range,  and  it  is  not  suitable,  e.g.  for  nonlinear 
spectroscopy. 

*In  conclusion  we  see  that  the  FEL  offers  the  possibility  of  a  source  with 
wide  continuous  tunability,  high  peak  power,  short  Fourier-transform-llmlted 
pulses,  that  are  superior  to  all  other  sources. 

IV.  Physical  Phenomena  in  the  FIR 

The  range  of  phenomena  of  interest  in  connection  with  the  FIR  spectrum  is 
very  wide.  For  example: 16 

— Superconductor  gaps  have  energies  corresponding  to  the  FIR. 

— "Hot"  electrons  in  GaAs/GaAlAs  junctions,  sub-bands  in  GaAs/GaAlAs 
superlattices  and  quantum  wells  are  phenomena  showing  FIR  resonances  and  fast 
decays  to  be  studied  with  short  pulses,  and  their  coherent  and  nonlinear 
behaviour  is  of  Interest  in  connection  with  possible  electronic  devices. 

— Plasmons  in  semiconductors  can  have  all  resonance  frequencies  from  the 
visible  to  the  FIR. 

— Low-energy  electronic  excitations  of  Impurities  in  crystal  lattice  show 
narrow  atomlike  lines. 

— Molecular  gas  rotations. 17 

— Intermolecular  vibrations  in  molecular  crystals. 

— Antiferromagnetic  resonances  are  often  in  the  sub-mm  range. 

V.  Possible  Uses  of  FIR  FEL 

Essentially,  the  whole  FIR  spectrum  can  be  covered  by  blackbody  sources  or 
C02-point  contact  diode  mixers.  But  there  are  essentially  two  aspects  that 
are  not  studied  by  these  sources,  and  only  in  a  limited  way  by  FIR  gas 
lasers: 

— High  (instantaneous)  intensities  (which  need  high  brightness)  of  FEL  allows 
excitation  of  physical  systems  in  the  nonlinear  regime. 

— High  intensity  coherent  (Glauber-state)  short  pulses  allow  coherent  excita¬ 
tion  of  physical  systems,  that  is  to  study  phenomena  that  can  not  be  described 
by  "rate  equations"  alone,  or,  in  terms  used  in  microwave  magnetic  resonance, 
to  measure  the  T2  (dephasing  time)  as  well  as  the  (energy  decay  time) . 

Then,  in  thinking  about  the  use  of  FIR  FEL,  one  should  not  only  think  to  the 
materials  that  have  FIR  resonances,  but  to  the  aspects  for  which  the  FEL  is 
unique;  nonlinearity  and  coherence. 
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It  is  probably  only  with  the  FEL  that  materials  can  be  studied  for  their 
possible  use  as  solid-state  FIR  sources,  mixers,  modulators  etc.,  and  this 
might  open  the  possibility  of  extending  the  microwave  electronics  to 
frequencies  of  the  order  of  a  terahertz  (1  THz  =  10l2  Hz) . 

In  the  last  two  years,  some  experiments  have  been  done  with  the  Santa  Barbara 
100-400  ym  FEL.^ 

VI.  Example:  the  Hefei  (USTC)  Linac  FEL 

At  the  USTC  in  Hefei  (Anhui) ,  an  experiment  is  under  way  to  study  a  linac- 
based  FEL  at  10  ym  wavelength.^®  At  present  the  linac  beam  characteristics 
have  still  to  be  improved  in  order  to  obtain  enough  gain  for  oscillation  at 
10  ym. 

It  is  interesting  to  see  how  operation  at  100  ym  would  give  a  higher  gain  and 
could  be  obtained  with  a  lower  beam  quality.  The  present  characteristics  are: 

Y  =  40-80.,  e  =  Sir  10~®  m.  rad,  Ay/y  -  10“^. 

With  an  undulator  of  98  periods  of  3  cm  and  K  =  0.8,  the  output  at  y  =  44  is 
X  =  10  ym. 

The  gain  is  very  low  because  eq.  4  is  not  satisfied  (and  some  reduction  comes 
also  from  eq.  5). 

At  100  ym  eq.  4  would  be  satisfied,  and  the  gain,  with  y  =  15,  would  be  more 
than  20  percent,  which  is  comfortable.  The  same  could  be  obtained  at  y  =  30 
with  an  undulator  with  N  =  75  and  Xo  =  10  cm. 

VII.  Conclusions 

It  seems  not  only  realistic  and  less  expensive  to  aim  first  at  realising  a 
FEL  in  the  FIR  (before  aiming  at  near  IR  or  UV) ,  but  also  there  is  a  particu¬ 
lar  interest  in  the  fact  that  many  nonlinear  and  coherent  phenomena  in  the 
FIR  can  be  studied  almost  only  with  the  FEL  (while,  e.g.  in  the  near  IR  and 
visible  there  are  high  peak  power,  tunable  lasers,  continuous -wave  or  short 
pulse  (Fourier-transform  limited),  as  well  as  mirrors,  modulators,  fibres,...), 
and  these  studies  could  be  important  for  the  possible  development  of  a  THz 
electronics . 
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